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Abstract 

We study the general properties of fluid spheres satisfying the 
heuristic assumption that theirs areal and proper radius are equal 
(the Euclidean condition). Dissipative and non-dissipative models are 
considered. In the latter case, all models are necessarily geodesic and 
a subclass of the Lemaitre-Tolman-Bondi solution is obtained. In the 
dissipative case solutions are non-geodesic and are characterized by 
the fact that all non-gravitational forces acting on any fluid element 
produces a radial three-acceleration independent on its inertial mass. 
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1 Introduction 



Analytical or numerical solutions to Einstein equations describing dissipative 
gravitational collapse are thought to be useful not only for describing specific 
astrophysical phenomena, but also as test-bed for probing cosmic censorship 
and hoop conjecture among other important issues. 

In this work we consider a large family of solutions, derived from the 
heuristic assumption that the areal radius of any shell of fluid, which is the 
radius obtained from its area, equals the proper radial distance from the cen- 
tre to the shell. Since these two quantities are always equal in the Euclidean 
geometry, systems described by solutions satisfying such a condition will be 
called "Euclidean stars". 

Some of the models are necessarily dissipative. This is appealing from 
a physical point of view, since gravitational collapse is a highly dissipative 
process (see [1, 2, 3] and references therein). This dissipation is required to 
account for the very large (negative) binding energy of the resulting compact 
object of the order of — 10 53 erg. 

The resulting dissipative models have a distinct dynamical property, namely 
all non-gravitational forces acting on any fluid element, produce a radial 
three-acceleration being independent of the inertial mass density of the fluid 
element. This behaviour, which is characteristic of the gravitational force, is 
now shared, due to the Euclidean condition, by all forces. The specific case 
of the shear-free and conformally flat fluid, are considered in detail. 

Non-dissipative models are necessarily geodesic, belonging to the Lemaitre- 
Tolman-Bondi (LTB) solutions (more specifically to the parabolic subclass). 
They may describe collapsing dust or, more generally, anisotropic fluids [4]. 

2 The Euclidean condition and its consequences 

We consider a spherically symmetric distribution of collapsing fluid, bounded 
by a spherical surface S. The fluid is assumed to be locally anisotropic with 
principal stresses unequal and undergoing dissipation in the form of heat 
flow. Choosing comoving coordinates inside S, the general interior metric 
can be written 

dsl = -A 2 dt 2 + B 2 dr 2 + R 2 (d9 2 + sin 2 Odtf), (1) 
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where A, B and R are functions of t and r and are assumed positive. We 
number the coordinates x° = t, x l = r, x 2 = 9 and a; 3 = <fi. 
The matter energy-momentum inside £ has the form 

T *p = (A* + P±)V a V p + P ± g af} + (P r - Pi_)x«Xp + q a Vp + V a qp, (2) 

where \i is the energy density, P r the radial pressure, P± the tangential 
pressure, q a the heat flux, V a the four-velocity of the fluid and x a a un it 
four-vector along the radial direction. These quantities satisfy 

V a V a = -l, V a q a = 0, X a Xa = l, X a V a = 0. (3) 
The four-acceleration a a and the expansion 9 of the fluid are given by 

a a = V a . p V p , e = V a ;a , (4) 

and its shear by 

V a p = V(a;/3) + d(aV/3) ~ ^9(^/3 + V a Vp). (5) 

Since we assumed the metric (1) comoving then 

V a = A~X, q a = qB- l 5«, x a = B- 1 6?, (6) 

where q is a function of t and r. From (4) with (6) we have for the four- 
acceleration and its scalar a, 

A' ( A' \ 2 

a 1 = — , a 2 = a a a a = — , (7) 



and for the expansion 



where the prime stands for r differentiation and the dot stands for differen- 
tiation with respect to t. With (6) we obtain for the shear (5) its non zero 
components 

^ D 2 °33 1 D 2 / n \ 

O-ll = o B ^22 = — T7 = "q a ' ( 9 ) 

3 sin 3 
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and its scalar 



where 



<? af3 <? a p = 3O- 2 



_ 1_ I B _ IV 
~ A\B~R, 



(10) 

(11) 



The mass function m(t, r) introduced by Misner and Sharp [5] (see also 
reads 



m 



R 3 

p 23 



R 
2" 



IV 



(12) 



We can define the velocity U of the collapsing fluid as the variation of the 
areal radius with respect to proper time, i.e. 



U = DtR < (in the case of collapse), 



(13) 



where D T = (l/A)(d/dt) defines the derivative with respect to proper time. 
Then (12) can be rewritten as 



E- R ' 



(14) 



The proper radial three-acceleration D T U of an infalling particle inside 
S can be calculated to obtain 



(in \ 
— + AixP r R \ + Ea, 



(15) 



feeding back this expression into the radial component of the Bianchi iden- 
tities produces (see [7] for details) 



-E 



(Tfl \ 
— + 4nP r R) 

D T q + 2q (J^ + a^j 



D R P r + 2(P r - P ± ) 



(16) 



The physical meaning of different terms in (16) is discussed in detail in [1, 
7, 8]. We would like just to recall that the first term on the right hand side 
describes the gravitational force term. As expected from the equivalence 
principle, its contribution to D T U is independent on the inertial mass density 
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[i + P r . The two last terms describe non-gravitational force terms (i.e. their 
combination vanishes in a geodesic motion). 

Two radii are determined for a collapsing spherical fluid distribution by 
the metric (1). The first is determined by R(t,r) representing the radius as 
measured by its spherical surface, hence called its areal radius. The second 
is obtained out its radial integration / B(t,r)dr, hence called proper radius. 
These two radii in general, in Einstein's theory, need not to be equal, unlike 
in Newton's theory. Here we assume those two radii to be equal. Hence with 
this condition we can write, 



implying from (14) 



B = R', 



E 



The field equations with this condition become 



kP\ = - 




(17) 
(18) 

(19) 

(20) 
(21) 



(22) 



while the mass function (12) now reads, 




(23) 



It is clear from (23) that if R = then m = and spacetime becomes 
Minkowskian. Therefore all Euclidean stars are necessarily non-static. Fur- 
thermore, using (13), (23) can be rewritten as 



m 
R 



2 ' 



(24) 
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Hence, (24) can be interpreted as the Newtonian kinetic energy (per unit 
mass) of the collapsing particles being equal to their Newtonian potential 
energy. 

From (20), we observe that if the system is dissipating in the form of heat 
flow, the collapsing source needs A' ^ 0, implying because of (7) a a ^ 0. This 
means that dissipation does not allow collapsing particles to follow geodesies. 
Inversely, of course, non-dissipative Euclidean models are necessarily geodesic, 
since q = implies because of (7) and (20) that a a = 0. 

It is interesting to observe that due to the Euclidean condition, the dy- 
namical equation (15) or (16) becomes, 



implying that the non-gravitational force term (the last on the right hand 
side) contributes to D t U, for any fluid element, independently on its in- 
ertial mass density. In other words, the Euclidean condition produces a 
"gravitational-like" behaviour in non-gravitational forces (which are con- 
trolled by q). Thus, the effect of non-gravitational forces amounts to modify 
the gravitational force term, leaving a "gravitational-like" force term pro- 
ducing a radial three-acceleration independent on the inertial mass density 
of the fluid element. 

The Weyl tensor C a p 7 s for metric (1) with (17) has the following non zero 
components, 




(25) 




(26) 



and all the other non zero components are proportional to (26), 




(27) 
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With (11), (21) and (22) we can rewrite (26) like 



c. 



AR 



12 



0101 



R 

k(P± -P r )A + 2—a 



(28) 



showing that for isotropic systems the shear-free conditions implies a confor- 
mally flat source. 

We consider next the non-dissipative case. 



3 Collapse with q = 

As mentioned before, for this case we have from (20) that A' = which 
means A — A(t) and by rescaling t we can have 

A—l. (29) 

Of course such models are members of the Lemaitre-Tolman-Bondi (LTB) 
spacetimes [9, 10, 11], furthermore they correspond to the parabolic case. 
Indeed, the general metric for LTB spacetimes read, 

ds 2 = -dt 2 + dr 2 + R 2 (d9 2 + sin 2 6d<p 2 ) , (30) 

1 — Kyr) 

where K{r) is an arbitrary function of r. 

Imposing the Euclidean condition (17) in (30), one obtains K — 0, which 
defines parabolic LTB spacetimes. Further, assuming that the source consists 
of pure dust (P r — P± — 0) then it follows from the field equations that 

R(t,r) = [ Cl (r)t + c 2 (r)f 3 , (31) 

and f 

Kfl = 3( Cl t + c 2 )(^ + c' 2 )' (32) 

where Ci(r) and c 2 (r) are integration functions. Hence the solution reduces 
to parabolic LTB collapsing dust [9, 10, 11]. 
From (29) and (31) we have for (26), 

m 4 Cl (c;t + c 2 ) 2 
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where a, from (11), is 



= c 2C ; - ci4 

(c 1 t + c 2 )(c' 1 t + c , 2 )' 1 ; 

In the shear- free case, C\ = c 2 , the system becomes conformally flat too, 
and with the freedom for choosing the r coordinate we can assume c\ = r 3//2 
recovering the Friedmann critical dust sphere. 

Of course more general models can be obtained by relaxing the condition 
of vanishing pressure, we recall that LTB spacetime is compatible with an 
anisotropic fluid [12, 13]. 



4 Collapse with q ^ 

We consider now the dissipative case. For simplicity we assume the fluid to 
be shear-free. In this latter case the line element can be written as [14] 



dsl = -A 2 dt 2 + B 2 [dr 2 + r 2 (d6 2 + sin 2 Od^ 2 )}, 
then the Euclidean condition becomes 

B = (Br)' ^B = fit) 

implying 

R = f(t)r, 

where / is an arbitrary function of t. 
The field equations (19-22) now read, 



l 

I 2 



4" 



A 
l A 



kP\ = - 



A 2 



4" 



A 



f f 



Kq - 

1)1 
f) f 

1_ fAT_ 

T 2 \A 



3 'i) 

f A' 



-2 



+ 



P A 2 '' 
2 A' 

r A I 



(35) 

(36) 
(37) 



(38) 
(39) 
(40) 
(41) 
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From (40) and (41) we have 

«(Pl -Pr) = ^ (f - . (42) 

From (28) it follows that (in the shear-free case) if the collapsing source is 
conformally flat, it must be isotropic in its pressures, and vice- versa. 

The general form of all conformally flat and shear-free metrics is known 
[15], it reads 



A 



ei (t)r 2 + l B, (43) 



where e\ is an arbitrary function of t, and 



B e 2 (t)r* + e 3 (ty {U) 

where e 2 and are arbitrary functions of t. 
The Euclidean condition then implies 

e 2 = 0, e 3 = j. (45) 

An approximate solution of this kind has been presented and discussed 
in [15]. Furthermore, an exact solution is also known [16], which in turn is a 
particular case of a family of solutions found in [17]. It reads (see Case III 
in [16]) 

f(t) = {fix + fofe- 2 ^ (46) 

and 

A = (ar 2 + 1)/, (47) 

where a, fix and f3 2 are constants. The above solution satisfies junction 
conditions and its physical properties have been discussed in [16]. Thus we 
shall not elaborate any further on it. Suffice to say at this point that its 
physical properties are reasonable and a thermodynamic analysis brings out 
the relevance of relaxational effects on the evolution of the system. 
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